TALBOT EFFECT FOR DISPERSION IN LINEAR OPTICAL FIBERS 

AND A WAVELET APPROACH 
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We shortly recall the mathematical and physical aspects of Talbot's self-imaging 
effect occurring in near-field diffraction. In the rational paraxial approximation, 
the Talbot images are formed at distances z = p/q, where p and q are coprimes, 
and are superpositions of q equally spaced images of the original binary trans- 
mission (Ronchi) grating. This interpretation offers the possibility to express the 
Talbot effect through Gauss sums. Here, we pay attention to the Talbot effect in 
q '■ the case of dispersion in optical fibers presenting our considerations based on the 

close relationships of the mathematical representations of diffraction and disper- 
q . sion. Although dispersion deals with continuous functions, such as gaussian and 

• i-H . supergaussian pulses, whereas in diffraction one frequently deals with discontin- 

uous functions, the mathematical correspondence enables one to characterize the 
Talbot effect in the two cases with minor differences. In addition, we apply, for 
the first time to our knowledge, the wavelet transform to the fractal Talbot effect 
in both diffraction and fiber dispersion. In the first case, the self similar character 
of the transverse paraxial field at irrational multiples of the Talbot distance is 
confirmed, whereas in the second case it is shown that the field is not self simi- 
lar for supergaussian pulses. Finally, a high-precision measurement of irrational 
distances employing the fractal index determined with the wavelet transform is 



> 

' pointed out 

o ' 
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1. Introduction 

Near field diffraction can produce images of periodic structures such as 



gratings without any other means. This is known since 1836 when this self- 
imaging phenomenon has been discovered by H.F. Talbot, one of the inventor 
pioneers of photography.^ Take for example a periodic object as simple as a 
Ronchi grating which is a set of black lines and equal clear spaces on a plate 
repeating with period a. In monochromatic light of wavelength A one can 
reproduce its image at a "focal" distance known as the Talbot distance given 
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by = o 2 A _1 as shown in Fig. 1. Actually this famous focal distance has 
been first derived by Lord Rayleigh in 18810 Moreover, more images show 
up at integer multiples of the Talbot distance. It was only in 1989 that the 
first and only one review on the Talbot effect has been written by Patorski.^ 
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In the framework of Helmholtz equation approach to the physical op- 
tics of a (5-comb grating, Berry and KleiiP showed in 1996 that the Talbot 
diffraction wavefield has a highly interesting arithmetic structure related to 
Gauss sums and other fundamental relationships in number theory. 

In this paper, after briefly reviewing the results of Berry & Klein, we 
show that analogous results can be obtained easily for the case of dispersion 
in linear optical fibers. Moreover, we apply for the first time in the literature 
the wavelet transform^! to the fractal Talbot problem. The point with the 
wavelet transform is that it contains more information with respect to the 
Fourier transform, which is behind the Helmholtz equation. Wavelet trans- 
forms have been considered as a very suitable mathematical microscope for 
fractals due to their ability to reveal the construction rules of fractals and 

to resolve local scaling properties as noticed before for the case of fractal 
aggregates ED 

2. Talbot effect for Helmholtz scalar diffraction fields 

The diffraction due to Ronchi gratings can be approached analytically us- 
ing the Helmholtz equation. Passing to dimensionless transverse and parax- 
ial variables £ = x/a and £ = z/a, respectively, the scalar wave solution 
C) of the Helmholtz equation can be expressed as a convolution in £ of 
the Ronchi unit cell square function 




»(0 - { o l j Vi in 

and the Dirac comb transmittance, i.e. 



/+1/2 / oo \ 

<?(£') E exp[i27m(£-£')]exp[ie n (C)] . 



(2) 



In the previous formulas, the unit cell is the single spatial period of the 
grating, which we take centered at the origin and of length equal to unity 



and n (C) = 27rC^-yT — (j^j is a phase produced by the diffraction of 
the Dirac comb 'diagonal' rays. The so-called Fresnel approximation for this 
phase, i.e., a Taylor expansion up to the second order for the square root, 
leads to & S (C) ~ —^n 2 ^. It can be easily shown now that in the Fresnel 
approximation Eq. |2]can be written as an infinite sum of phase exponentials 
in both variables £ and £ 

oo oo 

*p(£>0= E 3nexp[i27rn£ - mn 2 Q = ^ fn^p(CC), (3) 
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where the amplitudes g n are the Fourier modes of the transmittance function 
of the Ronchi grating 

/•+ 1 /4 , 

g n = d£ exp[-z2vrn£] . (4) 

J-l/4 

Furthermore, by discretizing the optical propagation axis £ by means of 
rational numbers, one can write the rational paraxial field as a shifted delta 
comb affected by phase factors, which is the main result of Berry and Klein: 



1 00 / n\ 
^72 6 Up )exp[iWn;«,p)] » ( 5 ) 

^ n=— 00 * ^ ' 

where £ p = £ — e p /2 and e p = 0(1) if p is even (odd). The paraxial phases 
exp[i<P ( j i fj r (n; q,p)] are specified in Section 4 and appear to be the physical 
quantities directly connected to number theory. At the same time, this 
rational approximation allows for the following physical interpretation of the 
paraxial self-imaging process: in each unit cell of the plane p/q, q images 
of the grating slits are reproduced with spacing a/q and intensity reduced by 
!/<?• 

3. Wave propagation in optical fibers 

Field dispersion phenomena in dielectric media are described in terms of 
wave equations with sources of the form 

1 d 2 E d 2 P 

-"w (6) 

that can be obtained from the Maxwell equations under minimal assump- 
tions on the constitutive equations. As known, in fiber optics technology, 
electromagnetic dispersion is defined in terms of the propagation constant 
(wavenumber) of each frequency mode @{uj) = n(u)^. In the following we 
will use one of the simplest optical fibers having a core-cladding step profile 
of the index of refraction. In addition, the famous slowly varying envelope 
approximation (henceforth SVEA) is a realistic approach when treating the 
propagation of quasi- monochromatic fields, such as laser fields and other 

types of coherent beams within such materials. For more details we refer 

171 

the reader to textbooks. ' 

SVEA means decomposing the electromagnetic fields in two factors: a 
rapidly varying phase component and a slowly varying amplitude field A 
enveloping the rapid oscillatory fields. The following Schrodinger-like dis- 
persion equation can be obtained for A in the SVEA approximation 
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where [3% is the second coefficient in the Taylor expansion of the propagation 
constant in the neighbourhood of the central resonance frequency. This 
is the simplest form of the dispersion equation that one can envision in 
which actually no material propagation shows up. It can be fulfilled in 
the practical situation when the dielectric medium has sharp resonances 
(5u r <C u> r ). Because of technological requirements, /?2 is usually a negative 
parameter corresponding to the so-called anomalous dispersion region. As 
can be seen, the SVEA equation has exactly the same mathematical form as 
the diffraction equation in the paraxial approximation: 

where is the electric field close to the propagation axis. 

4. From diffraction to fiber dispersion 

Many results in diffraction can be translated to the case of dispersion in 
fibers by using the following substitutions 

x — > t 
y -> r 
z — > z. 



In the first row one passes from the grating axis to the time axis of a frame 
traveling at the group velocity of a pulse. In the second row one passes from 
the second grating axis that here we consider constant to the transverse 
section of the optical fiber. Finally, the propagation axis remains the same 
for the two settings. This change of variables will be used here to compare 
the results obtained in the two frameworks. 

The general solution of the SVEA dispersion equation (JJJ) for the am- 
plitude A(z, t) depends on the initial conditions. Assuming a periodic input 
signal written as a Fourier series, i.e., A(0,t) = I]n=^°oo C°e~ w ™*, where C° 
are the Fourier coefficients of the initial pulse at the entrance of an opti- 
cal fiber with linear response, one can write the pulse at an arbitrary z as 
follows: 

where u n = 2im/T. (9) 

If the scaling of variables r = t/T, ( = 2z/zt is employed, A(z, t) can be 
rewritten as 

A(C,r) = ^C°exp 



A(z,t) = ^C°exp 



iuj n t 



mn £ — i2nnT 



(10) 
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because the Talbot distance corresponding to this case is zt = T 2 /ir. Just 
as in the context of diffraction, the convolution of the unitary cell with the 
propagator can be equally done before or after the paraxial approximation 
is employed. We notice that Eq. can be also written as 



MC,r) 



T/2 
-T/2 



A(0,T')a((,T f -T)dr' 



(11) 



since C° are nothing but the Fourier coefficients of the input signal and 
where 



q(C,t)= ^ exp iim 2 £ — i2nnT 



(12) 



can be thought of as the analog of the paraxial propagator El In this expres- 
sion, the trick is to turn the continuous propagation axis into the rational 
number axis and also to perform the integer modulo division of n with re- 
spect to the rational denominator of the propagation axis, i.e., 



c 



p 
q 



n 



Iq + s. 



(13) 



Through this approximation, the sum over n is divided into two sums: one 
over negative and positive integers I, and the other one over s = n (modq) 



a(C,r) 



oo 3—1 
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iir(lq + s) 2 i2ir(lq + s)i 



(14) 



This form of r) is almost exactly the same as given by Berry & KleiiP 
and by Matsutani and Onishi0The difference is that the sign of the exponent 
is opposite. Following these authors one can express a in terms of the Poisson 
formula leading to 



2 oo ' ^ (J— i 

Q(< ' T)= vf„EbS exp 
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S(r p + -) , (15) 

q 



where t p is a notation similar to £ p . The rest of the calculations are straight- 
forwardly performed though they are lengthy. By algebraic manipulations 
the phase factor can be easily obtained and we reproduce below the two 
expressions for direct comparison 



$di sp (n;q,p) 



1 9_1 f 7T 



ps 2 — 2s(n 



qe P , 



(16) 
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diflfr 



1 



q-1 



n;q,p) = ^ exp 



s=0 



TV 
I — 

q 



ps 



(17) 



Both phases are special types of Gauss sums from the mathematical stand- 
point (see the Appendix). The difference of signs here appears because of 
the sign convention chosen for the Fourier transform. Not surprisingly, the 
changes in the mathematical formulation are minimal although the experi- 
mental setup is quite different. 

If one tries to make computer simulations using the Fourier transform 
method, the Gibbs phenomenon is unavoidable for discontinuous transmit- 
tance functions. However, in the case of fiber dispersion, one class of contin- 
uous pulses one could work with are the supergaussian ones, i.e., functions 
of the following form 



A(C = 0, r) = A exp 



_v 



O"0 



(18) 



where TV" is any even number bigger than two. The bigger the chosen N 
the more the supergaussian pulse resembles a square pulse. A computer 
simulation of the evolution of a supergaussian pulse train is given in Fig. 2.^ 



5. Irrational Talbot distances 



5.1. Fractal approach 

In the Talbot terminology the self-reconstructed images in the planes 
z = (p/q)zT consist of q superposed copies of the grating as already men- 
tioned, completed with discontinuities. Although there is a finite number of 
images at fractional distances, they still represent an infinitesimal subset of 
all possible images that occur at the irrational distances. 

However, in the planes located at irrational fractions of the Talbot dis- 
tance the light intensity is a fractal function of the transverse variable. The 
field intensity has a definite value at every point, but its derivative has no 
definite value. Such fractal functions are described by a fractal dimension, 
D, between one (corresponding to a smooth curve) and two (corresponding 
to a curve so irregular that it occupies a finite area). In the case of Ronchi 
gratings, for example, the fractal dimension of the diffracted intensity in the 
irrrational transverse planes is 3/2.^^ 

Moreover, one can define the so-called carpets which are wave intensity 
patterns forming surfaces on the (£ , £) plane, i.e., the plane of the transverse 
periodic coordinate and the propagation coordinate. Since they are surfaces, 
their fractal dimension takes values between two and three. According to 
Berry, in general for a surface where all directions are equivalent, the fractal 
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dimension of the surface is one unit greater than the dimension of an inclined 
curve that cuts through it. Taking into account that for a Talbot carpet the 
transverse image curves have fractal dimension D = 3/2, then the carpet's 
dimension is expected to be 5/2. However, Talbot landscapes were found not 
to be isotropic. For fixed £, as the intensity varies as a function of distance £ 
from the grating, the fractal dimension is found to be 7/4, one quarter more 
than in the transverse case. Therefore the longitudinal fractals are slightly 
more irregular than the transverse ones. In addition, the intensity is more 
regular along the bisectrix canal because of the cancelation of large Fourier 
components that has fractal dimension of 5/4. The landscape is dominated 
by the largest of these fractal dimensions (the longitudinal one) , and so is a 
surface of fractal dimension 1+ 7/4 = 11/4. 

5.2. Wavelet approach 

Wavelet transforms (WT) are known to have various advantages over the 
Fourier transform and in particular they can add up supplementary infor- 
mation on the fractal features^ A simple definition of the one-dimensional 
wavelet transform is 

/oo 
/(s)/4 )n (s)ds (19) 
-oo 

where h mn is a so-called daughter wavelet that is derived from the mother 
wavelet h(s) by dilation and shift operations quantified in terms of the dila- 
tion (m) and shift (n) parameters: 

1 , / s — n\ . . 

/l m ,n S = -=h 20 

ym V rn / 

In the following, we shall use the Morlet wavelets which are derived from 
the typical Gaussian-enveloped mother wavelet which is itself a windowed 
Fourier transform 

h(s) = exp[-(s/s ) 2 ] exp(z27rfcs). (21) 

The point is that if the mother wavelet contains a harmonic structure, e.g., 
in the Morlet case the phase exp(i2irks)), the WT represents both frequency 
and spatial information of the signal. 

In the wavelet framework one can write the expansion of an arbitrary 
signal (p(t) in an orthonormal wavelet basis in the form 

^)=EE«) . (22) 

m n 
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i.e., as an expansion in the dilation and translation indices, and the coeffi- 
cients of the expansion are given by 

/oo 
tp{t)W mjn (t)dt . (23) 
-oo 

The orthonormal wavelet basis functions W m n(t) fulfill the following dilation- 
translation property 

W m , n {t) = 2 m / 2 W(2 m t - n) . (24) 

In the wavelet approach the fractal character of a certain signal can be 
inferred from the behavior of its power spectrum P(f), which is the Fourier 
transform of the autocovariance (also termed autocorrelation) function and 
in differential form P{f)df represents the contribution to the variance of a 
signal from frequencies between / and / + df. Indeed, it is known that for 
self-similar random processes the spectral behavior of the power spectrum is 
given b y 11 

P^uj) ~ \uj\~^' , (25) 

where 7/ is the spectral parameter of the wave signal. In addition, the 
variance of the wavelet coefficients possesses the following behavioJ^ 

var<«(2 m )^/ . (26) 

These formulas are certainly suitable for the Talbot transverse fractals 
because of the interpretation in terms of the regular superposition of identical 
and equally spaced grating images. We have used these wavelet formulas in 
our calculations related to the same rational paraxiallity for the two cases 
of transverse diffraction fields (Fig. 3) and the fiber-dispersed optical fields 
(Fig. 4), respectively. The basic idea is that the above-mentioned formulas 
can be employed as a checking test of the self-similarity structure of the 
optical fields. The requirement is to have a constant spectral parameter 
7j over many scales. In the case of supergaussian pulses, their dispersed 
fields turned out not to have the self-similarity property as can be seen by 
examining Fig. 4 where one can see that the constant slope is not maintained 
over all scales. In Figs. 5 and 6 the behavior of the wavelet transform 
using Morlet wavelets for the diffraction field is displayed. A great deal 
of details can be seen in all basic quantities of the diffracted field, namely 
in the intensity, modulus, and phase. On the other hand, the same wavelet 
transform applied to the N=12 supergaussian dispersed pulse (see Fig. 7), 
although showing a certain similarity to the previous discontinuous case, 
contains less structure and thus looks more regular. This points to the fact 
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that if in diffraction experiments one uses continuous transmittance gratings 
the fractal behavior would turn milder. 

More realistic ally , paraxial waves display electric and magnetic polariza- 
tion singularities!*^ If the paraxial wavefield is treated as a signal then it is 
worth pointing out here that detection of signal singularities has been stud- 
ied in quite detail by the experts in wavelet processing!^ We plan to study 
this aspect in future research. 

6. Conclusion 

The fractal aspects of the paraxial wavefield have been probed here by 
means of the wavelet transform for the cases of diffraction and fiber disper- 
sion. In the case of diffraction, the previous results of Berry and Klein are 
confirmed showing that the wavelet approach can be an equivalent and more 
informative tool. The same procedure applied to the case of fiber dispersion 
affecting the paraxial evolution of supergaussian pulses indicates that the 
self-similar fractal character does not show up in the latter type of axial 
propagation. This is a consequence of the continuous transmittance func- 
tion of the supergaussian pulses as opposed to the singular one in the case 
of Ronchi gratings. 

Finally, as a promising perspective, we would like to suggest the following 
experiment by which irrational distances can be determined. The idea is 
that the spectral index of the Talbot fractal images can be used as a very 
precise pointer of rational and irrational distances with respect to the Talbot 
one. Suppose that behind a Ronchi grating under plane wave illumination a 
CCD camera is mounted axially by means of a precision screw. The Talbot 
image at can be focused experimentally and can be used to calibrate the 
whole system. An implemented real time wavelet computer software can 
perform a rapid determination of the fractal index jf, which in turn allows 
the detection of changes of the distance in order to determine if the CCD 
camera is at rational or irrational multiples of the Talbot distance. To the 
best of our knowledge, we are not aware of another experimental setup in 
which irrational distances can be determined in such an accurate way. This 
also points to high-precision applications in metrology. 
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Appendix: Gauss sums 

It is practically trivial to show that 

6-1 

£ e^i = 0, 

n=0 

for any integer b bigger than one. However, similar sums of nonlinear powers 
in n are usually nonzero 

b- 1 k 

£ e°*V / 0, 

n=0 

where > 2, b is a prime with 6=1 (mod fc) and the sum is over an 
arbitrary complete system of residues (mod b) . Such sums are really difficult 
to calculate. For k = 2, they are known as (quadratic) Gauss sums because 
two centuries ago Gauss first gave explicit results for them but only after 
4 years "with all efforts" . A paper by Berndt and EvanJ^ is an excellent 
survey of the history of Gauss sums, see also the appendix of the work of 
Merkel et al in this volume. Gauss sums having linear phases added to 
the quadratic ones are of special interest in physical applications such as 
the optical ones discussed in this paper. We concentrate here on this type 
of Gauss sums following the papers of Hannay and Berry-L^ anc [ a i so of 
Matsutani and Onishi)^ 

At rational values of the propagation coordinate the diffraction and dis- 
persion phases can be expressed through the infinite Gaussian sum 



I Nab 



G 00 (a,b,c) = Jim - uat > ex P 



N^oo 2abN v . 

m=~N ab 



'l(< 



i— ( am 2 + cm 



which is an average over a period 2abN. If a and b are co-prime numbers, 
the sum is zero unless c is an integer. Let m = bn + s, then the sum is 
divided into two sums 



Goo(a,6,c) = lim — — V V exp i\ (a{bn + s) 2 + c(bn + s] 

n=—Na s=U 

The argument of the exponential in the sum becomes 



1 b 



— (a(bn + s) 2 + c(bn + s)\ = ilixans + mn{abn + c) + i— (as 2 + bs 



and if ab and c are both odd or even at the same time, only terms in s 
survive: if n is even, the term is even and if n is odd, i.e. n = 2k\ + 1, and 
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also ab = 2k 2 + 1 and c = 2k 3 + 1, the term is 

mn{abn + c) = iw(2ki + l)[{2k 2 + l)(2fci + 1) + (2k 3 + 1)] 
= mt(2/ci + 1)(4£;i£;2 + 2/c2 + 2/ci + 2/c3 + 2) 
= i2vr(2/c 1 + l)(2/c 1 A;2 + /c2 + fei + /c3 + l) 



Thus Goo is rewriten as 



Na 6-1 



Goo (a, c) = lim 



ex p *t f a* 2 + 

ftTV ^ ft v 



r.s- 



n=-JVas=0 

and taking the limit, Goo coincides with the quadratic Gauss sum 

1 a ^ r 7r 

Goo(a, c) = S(a, ft, c) = — ^ exp i— |am 2 + cmj 

rrt=0 

Notice the following property: in the normalization fraction and the upper 
limit of the sum, a can be interchanged for b simultaneously. We will use 
this property because it is appropriate from the standpoint of the physical 
application to run the sum over the denominator of the phase, i.e., over ft. 

We consider now the various cases according to the odd even character 
of the parameters a, 6 and c. Only the first case will be treated in detail. 

Case a and c even, ft odd (because of coprimality to a) 

A factor exp (mcmk) can be added since c is even. Then 

1 6-1 

S(a,b,c) = -J2 e * fam2 e i T cm ( kb + 1 ) 



m=0 



Noticing that kb+1 = 1 (mod ft) and defining aa& = 1 (mod ft) we can write 

^ 6-i /[ a 1\ 
S(a, ft, c) = - ^2 ex P f ^r- ( m2 + cmabj J 

rn=0 ^ ^ J ' 

and completing the square in the argument of the exponential: 



6-1 



S(a,b,c) = - exp 



ITT- 



m=0 



m + 



ca&Y 
2 J 



cat, 
~2 ' ) 
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Next, one should notice the following simpler form: 

21 
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ex p 

m=0 



ivr— m + 



2 ) 



6-1 
^ exp 



?1=0 



m—n' 
ft 
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since the linear term plus the independent one in the left hand side expo- 
nential lead to a different order of the terms with respect to the right hand 
side of the equation. Furthermore we note that 



6-1 



5>xp 



n=0 



a n 
nr—n 
b 



6-1 

£ 

n=0 



1 + 



an/2 



exp ( nr—n 



where denotes the Legendre symbol that for 1% an odd prime number 

is defined as zero if l\ is divided by I2, 1 if h is a quadratic residue mod I2 
(there exists an integer k such that k 2 = l\ (mod I2)), and -1 if l\ is not a 
quadratic residue. 

Therefore the Gauss sum is expressed as follows: 



S(a, b, c) 



exp ( -iir-(ca b /2) 2 ) ^ 

n=0 



6-1 



1 + 



an/2 



exp [ nr—n 



and the latter sum under splitting into two sums yields zero in the first 
one. Making use of the Legendre symbol multiplicative property in the top 
argument we get: 



1 f a/2\ ( a 2 \ { n \ ( n 

S{a,b,c) = - ^— j ^exp \^-m-(ca b /2) j ^ ^-J^exp \^2ir- 



The following result from number theory*^ should be used in order to 
get the most compact result for S(a,b,c): 

l/2\ (2\ , ,,^-1)78 „_fb 2 -l\ 



and then 

S© i ex K^) = ^ exp (^^ ! ) 

for b = 1 (mod 4) or b = 3 (mod 4). Then: 

1 f"\ I , . e 2 - 1) . a, ,„., . (6 - 1) 2 \ 
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The terms in the exponential with the common factor i7r/8 are rewritten 

z2vr 



ITT 



± b 2 Tl + b i-2b + l)=- — (b-l) 
' 8 



having chosen the lower sign (there is no final change if the + sign is picked 
up). 

Recalling now that a and c are even we get 



S(a, b, c) 



1 fa 



Vb \b 
Case a, b, and c odd 



exp —ITT 



i(6-l) + ^( C a 6 /2) 2 



Exactly the same expression can be worked out with the only minor modi- 
fication that a is replaced by 2%a. 
Case a odd, b even, c odd 

Similar arguments lead to a slightly different result. 
The three final results are as follows: 



S(a, b, c) 



1 fa 
Vb \b 



exp 



(6-1) + -(ca b ) 2 
b 



a even, b odd, c even, 



S(a, b, c) 



1 fa 



exp 



.7T 



2a -, 



(6-1) + —% (ca b ) 



a odd, b odd, c odd, 



S*(a, b, c) 



-a + ^(ca b ) 



a odd, b even, c odd. 



These results were obtained under the assumption that b is prime. For 
the general case where b is not a prime the Legendre symbol should be 
substituted by the Jacobi symbol, which is a product of Legendre symbols 
defined in the following way. Let the prime decomposition of b be b = 
YliZiPi 1 - Then, by definition, the Jacobi symbol is 



n 



a 

Pi/ L 



In this way, the evaluation of the Talbot phases are performed by identi- 
fication of the parameters: a = ±p, b = q, c = +2(n =p qe p /2), respectively, 
where e p = (1) if p even (odd). The upper sign is for dispersion and the 
lower for diffraction. Thus, we get: 
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p even, q odd: 



$dis P (™;P, q) 
$diSr(n;p, q) 

p odd, q odd: 
®dis P (n; p,q) = 
$diSr(n;p,q) = 

p odd, q even: 
$dis P (™; P,q) 
®diffr(n; p,q) 



/ xp ri 
p\ /■ , 

exp +z— 



P 

£) 



cxp 

exp I -M 



7T 

'I 

7T 

4 



(g-l) + ^(2np,) 

q 

( g -l) + E(2np,) 
(g-l)+22^((2n-g)p 9 ) s 



(g-l) + 22^((2n + g)p^ 













' J 6xp 




— p + 


p> 














— p + 




l^exp 




p> 









-p+ -((2n + q)pg 
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Figure captions 



FIG. 1: Image of a Ronchi grating as obtained on a photographic plate 
located at the Talbot distance zt = 28.4 cm. This and many similar images 
have been obtained in the graduate work of Treviho- Gutierrez-^ in which a 
He-Ne laser working at its usual wavelength A = 632.8 nm has been used. 

FIG. 2: Computer simulation of |^4(C,t)| 2 of a N = 12 , do = 1.5 super- 
gaussian pulse train as given in Eq. (|18|) in a linear fiber characterized by 
the Schrodinger-like dispersion relation given in Eq. (JJJ). 

FIG. 3: (a) The fractal Talbot light intensity |^ p | 2 at ( = 144/377 and (b) 
the plot of the logarithmic variance of its wavelet coefficients (Eq. (|23|)). The 
line of negative slope of the latter semilog plot indicates fractal behaviour of 
the diffraction wavefield as we expected. The fractal coefficient is given by 
the slope and its calculated value is 7/. 

FIG. 4: Snapshot of the dispersed supergaussian pulse for N = 12 at 
C = 144/377 (close to the Golden Mean ( G = 3/2 - y/h/2). The log variance 
plot is monotonically decreasing displaying a plateau indicating a nonfractal 
behaviour of the N = 12 supergaussian pulse train. 

FIG. 5: The wavelet transform of the intensity |^ p | 2 at £ = 144/377 for 
(a) the unit cell and (b) half-period displaced grating unit cell. There is no 
difference because the square modulus is plotted. 

FIG. 6: Wavelet representations of: (a) the squared modulus of the ampli- 
tude and (b) phase of the Talbot diffraction field for fixed Q = 144/377 and 
a displaced unit cell. 

FIG. 7: Wavelet representations of the (a) amplitude and (b) phase of the 
Talbot dispersed supergaussian field [N = 12) for Q = 144/377. 
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